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We extend the renormalization of the NN interaction with Chiral Two Pion Exchange Potential 
to the calculation of non-central partial wave phase shifts with total angular momentum j < 5. 
The short distance singularity structure of the potential as well as the requirement of orthogonality 
conditions on the wave functions determines exactly the number of undetermined parameters after 
renormalization. 
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I. INTRODUCTION 



The original proposal by Weinberg 0, |3| and carried 
out for the first time by Ray, Ordoiiez and van Kolck fall, 
of making model independent predictions for NN scatter- 
ing using Chiral Perturbation Theory (ChPT) has been 
followed by a wealth of works 0, la, la, H IM IS flOl fllL [l2l 
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(for a review see e.g. Ref. 29]). The renormalized po- 
tential as given in Refs. 0,i4;i5j in configuration space is 
expanded taking m^/167r^/^ and m/M as small param- 
eters (m and M are the pion and nucleon masses respec- 
tively and / is the pion weak decay constant), with mr 
fixed. In this counting for the potential and in a given 
partial wave (coupled) channel with good total angular 
momentum the reduced potential can schematically be 
written as 



U{r) = Mm 



TO* TO 



M/(3)(r 



(1) 



where are known dimensionless functions which arc 
everywhere finite except for the origin and depend on 
the axial coupling constant. W^^^ depends also on three 
additional low energy constants ci — ciM, C3 — c^M 
and C4 = C4M which have been determined from ttN 
scattering ChPT studies in a number of works jsH Isil 

mm. 

At the level of approximation of Eq. Q these poten- 
tials are local and energy independent and become sin- 
gular at the origin. Thus, non-perturbative renormaliza- 
tion methods must be applied to give a precise mean- 
ing to the scattering amplitude ^341 (for a comprehensive 
review in the one channel case see e.g. Ref. [S^ and 
Ref. liSl for a modern perspective). Several methods 
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have be en prop osed to study the LO term in Eq. 0J for 
central l87|. [s^ lioL l4ll | and non-central waves. 
Recently |4ll |43 | we have shown how a renormalization 
program can be carried out for the NN interaction for 
the One Pion Exchange (OPE) and chiral Two Pion Ex- 
change (TPE) potentials in the central ^Sq and ^Si —^Di 
waves and its implications for the deutcron and pion- 
deuteron scattering |44j . In the present work we extend 
our analysis to all remaining partial waves with j < 5 
both for the OPE as well as for the chiral TPE poten- 
tials. As we showed in Refs. [H ^| the short distance 
behaviour of the chiral NN potential, Eq. , determines 
exactly how many counterterms are needed in order to 
generate renormalized and finite, i.e. cut-off indepen- 
dent, phase shifts. These counterterms can be deter- 
mined by fixing some low energy parameters while the 
cut-off is removed. It has been assumed that dimen- 
sional power counting in the counterterms can be made 
independently on the short distance singularity of the po- 
tential. This yields conflicts between naive dimensional 
power counting and renormalization which have been re- 
ported recently even for low partial waves |Q . So, one is 
led to an alternative: either one keeps the power count- 
ing and a finite cut-off or one removes the cut-off at the 
expense of modifying the power counting of the short 
distance interaction. The finite cut-off route has been 
explored in great detail in the past H H S H H ITol| . 
In this paper we explore further the possibility of tak- 
ing the alternative suggested by renormalization and the 
tight constraints imposed by finiteness. The analysis 
becomes rather transparent in coordinate space, where 
the counterterms can be mapped into boundary condi- 
tions m li^ at the origin. In practice renormaliza- 
tion may be carried out in several ways. In coordinate 
space it seems natural to exploit the locality of the long 
distance (renormalized) potentials and then to renormal- 
ize the full scattering problem. In the present work we 
adhere to this two-steps renormalization which has the 
additional advantage of making possible to determine a 
priori and based on simple analytical arguments the ex- 
istence of the renormalized limit and how many indepen- 
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dent renormalization conditions (counterterms) are com- 
patible with this limit. In this regard, let us remind that 
the main advantage of renormalization is that identical 
finite and unique results should be obtained regardless 
of the method of calculation (coordinate or momentum 
space) and regularization provided the same input phys- 
ical data are used to eliminate the divergencies. In par- 
ticular we also expect independence on the way how the 
limit is taken. 

The origin of the conflict can be traced back to the 
question whether for a given energy independent local 
potential, such as Eq. ^ , one can assume any short dis- 
tance physics regardless on the form of the long range 
potential. Renormalization group invariance, however, 
requires that any physical parameter sits on a renormal- 
ization trajectory and the corresponding evolution on the 
renormalization scale is dictated by the form of the long 
distance potential at all distances. The precise trajec- 
tory is uniquely fixed by a renormalization condition at 
very long distances. Thus, the separation between the 
short and long distance contribution is not only scale de- 
pendent but also potential dependent [s^ ^3- Renor- 
malization conditions are physical and do not exhibit 
this dependence. Finiteness of the scattering amplitude 
and orthogonality of scattering (and eventually bound 
state) wave functions impose very tight constraints on 
the allowed number of counterterms and their possi- 
ble scale dependence l^j ■ The discussion becomes 
rather straightforward in coordinate space and in terms 
of boundary conditions for ordinary differential equa- 
tions. In addition, unlike momentum space treatments, 
a very natural hierarchy of the renormalization problem 
takes place in configuration space |^ ^3 • More specifi- 
cally, orthogonality of different energy solutions requires 
an energy independent boundary condition on the wave 
function for the long distance local and energy indepen- 
dent potentials as it is the case for Eq. valid to NNLO, 
so that in all cases the effective range, and higher order 
threshold parameters cannot be taken as independent in- 
put parameters ^. 

The results found in Refs. 0, ^ can be concisely 
summarized as follows in the one channel case. For a 
regular potential, i.e., diverging less strongly than the 
inverse square potential, r^|f7(r)| < oo, one may choose 
between the regular and irregular solution. In the first 
case the scattering length is predicted while in the sec- 
ond case the scattering length becomes an input of the 
calculation. Singular potentials at the origin, i.e. ful- 
filling, r^|[/(r)| oo, do not allow this choice. If the 
potential is repulsive, the scattering length depends on 
the potential while for an attractive potential the scatter- 
ing length must be chosen as an independent parameter. 



^ Actually, the potential in Eq. contains distributional contri- 
butions, which strictly speaking are zer o for any finite distance. 
See the discussion in our previous work l43l 



In the coupled channel situation one must look at the 
strongest singularity of the potential eigenvalues at the 
origin, and apply the single channel results. 

In our formulation of the NN renormalization prob- 
lem threshold parameters play an essential role. Unfortu- 
nately, scattering threshold parameters for higher partial 
waves other than the S-waves have never been consid- 
ered in the context of chiral potentials H, IE S IE ^| ■ 
Instead, some calculations adjust their counterterms to 
fit the phase shifts in the region above threshold to the 
Nijmegen database ^3, ^ recent work we have 

filled the gap by carrying out a complete determination 
of these threshold parameters for the Reid93 and Nijmll 
potentials ■ the light of this new information it is 
quite possible that the good fits in the intermediate en- 
ergy region imply a somewhat less accurate description 
in the threshold region. This issue will become relevant 
in the description of some partial waves. 

The paper is organized as follows. In Sect.^we review 
the formalism for coupled channel scattering in the pres- 
ence of singular potentials at the origin. For completeness 
we list the potentials in Appendix IXI Based on the short 
distance behaviour of those potentials (see Appendix O 
and the requirement of orthogonality we determine the 
number of independent parameters for any partial wave 
with j < 5. In Sect. IIIII we present our results for the 
phase shifts. Specifically, we make a thorough analysis of 
cut-off dependence in all partial waves both for the OPE 
as well as for the chiral TPE potential. We also discuss 
the perturbative nature of peripheral waves within the 
present non-perturbative approach. Finally, in Sect. IIVI 
we present our conclusions. 



II. FORMALISM 

We solve the coupled channel Schrodinger equation for 
the relative motion which in compact notation reads. 



u"(r) + 



U(r) + ^ 



u(r) = fc^u(r) , 



(2) 



where U(r) = 2^„pV(r) is the coupled channel matrix 
reduced potential with /i„p = MpMn/{Mp + M„) the 
reduced proton-neutron mass which for j > can be 
written as. 



U°J(r) 
Ui^'(r) 



oj 



33 



fu'U,,-iir) U. 




(3) 







In Eq. lO F = diag(Zi(;i + 1),...,In{In + 1)) is the 
orbital angular momentum, u(r) is the reduced matrix 
wave function, k the CM. momentum and j the total 
angular momentum. In our case = 1 for the spin 
singlet channel with I — j and N — 3 for the spin triplet 
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channel with = j — 1, I2 = j and Z3 = j + 1. The 
otentials used in this paper were obtained in Refs. 0,0, 
in coordinate space and are hsted in Appendix IXI for 
completeness. 



In the present paper zero energy scattering parameters 
play an essential role since they are often used (see below) 
as input parameters of the calculation of phase shifts. 
Due to unitarity of the S- matrix in the low energy limit, 
k —^ 0, we have 



A. Long distance behaviour 

At long distances, we assume the usual asymptotic nor- 
malization condition 



u(r) -> h(-)(r) -h(+)(r)S, 



(4) 



with S the coupled channel unitary S-matrix. The corre- 
sponding out-going and in-going free spherical waves are 
given by 



h(±)(r) - diag(/i±(fcr),...,/i± (fcr)), 



(5) 



with h^{x) the reduced Hankel functions of order /, 



hi {x) = xH^^^^^{x) ( {x) 



), and satisfy the 



free Schrodinger's equation for a free particle. 

For the spin singlet state, s = 0, one has I — j and 
hence the state is uncoupled 



(6) 



whereas for the spin triplet state s = 1, one has the 
uncoupled / — j state 



(S - 1);, ; = -2ia;,,;fc' 



I'+i+i 



(11) 



with am the (hermitian) scattering length matrix ^. The 
threshold behaviour acquires its simplest form in the 
SYM representation. 



(12) 

(13) 
(14) 
(15) 
(16) 



In the BB form one has similar behaviours for the (5's but 
for Co which behaves as fc^^ instead of k'^^'^^ 



5^' 



1 



(17) 
(18) 

(19) 



S 



(7) 



B. Short distance behaviour 



and the two channel coupled 1,1' = j ± 1 states for 
which we use Stapp-Ypsilantis-Metropolis (SYM or Nu- 
clear bar) 49j parameterization 



( cos(2ej)e2*'^)-i i sin (2ej)e*(*^'-iH'+i) \ 
\ isin(2ej-)e'('^)-i+'^)+i) cos (2ej-)e2«5jii j 

In the discussion of low energy properties we also use 
the Blatt-Biedenharn (BB or Eigen phase) parameteri- 
zation defined by 



S^3 



cos 

sin e, 



cos e," 



— sme 
cose 



sme 



sme, cose 






e^'^iii ^ 



(8) 



The relation between the BB and SYM phase shifts is 

(9) 
(10) 



sin ((5 



i-i 



tan(2ej) 



tan(2ej) 



The form of the wave functions at the origin is uniquely 
determined by the form of the potential at short distances 
(see e.g. 0,IH| for the case of one channel and 0,|43| for 
coupled channels). For the chiral NN potential, Eq. Q), 
one has 



MC 



3,LO 



MC 



5,NLO 



MC 



6,NNLO 



ULo(r) 
Unlo('') 

UNNLo(f') 

(20) 

where LO includes the first term in Eq. NLO the 
first two terms and so on. Note that higher order po- 
tentials become increasingly singular at the origin. For a 
potential diverging at the origin as an inverse power law 



U(r) 



MCr 



(21) 



For non S-wave scattering the dimension of O; ;/ is fm'+' '^^ 
which is not a length. For simplicity we will abuse language 
and call them scattering lengths. 
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TABLE I: Sets of chiral coefficients considered in this work. 



Set 


Source 


ci(GeV"-') 


C3(GeV-^) 


C4(GeV"-') 


Set I 


TvN [31j 


-0.81 


-4.69 


3.40 


Set II 


NN [5] 


-0.76 


-5.08 


4.70 


Set III 


NN [7] 


-0.81 


-3.40 


3.40 


Set IV 


NN [8] 


-0.81 


-3.20 


5.40 



with C„ a matrix of generalized van der Waals coeffi- 
cients and n > 2 an integer. One diagonalizes the matrix 
C„ by a constant unitary transformation, G, yielding 



A/C„ Gdiag(±i?'i' 



,±R 



G 



(22) 



with Ri constants with length dimension. The plus sign 
corresponds to the case with a positive eigenvalue (re- 
pulsive) and the minus sign to the case of a negative 
eigenvalue (attractive). Then, at short distances one has 
the solutions 



u(r) ^ G 



'"i,±('')' 



(23) 



where for the attractive and repulsive cases one has 

Tl/4 



a. 



sm 



n~2 \ r 



+ (fi 



i/4 



exp 



Ri 
2 [ r 



(24) 
(25) 



respectively. This behaviour of the wave functions near 
the origin is valid regardless of the energy, provided the 
distances are small enough Here, (pi are arbitrary short 
distance phases which in general depend on the energy. 
There are as many short distance phases as short dis- 
tance attractive eigenpotentials. Orthogonality of the 
wave functions at the origin yield the relation 



N 
i=l 



r=0 



A 

E 



cos{ipi{k) - (pi{p)) , 



(26) 



where yl < iV is the number of the short distance attrac- 
tive eigenpotentials. 

The simplest choice to fix relative phases for a positive 
energy scattering state is to take the zero energy state 
p = as a reference state, and the zero energy short 
distance phase. In the particular case where only one 



^ In fact, the next correction to the near-the-origin wave func- 
tions, which is energy dependent, is suppressed by a relative 
(kR)'^ (r / R)"/"^^^ power with respect to the main term, so it is 
negligible in the r — > limit. 



eigenvalue is negative the short distance phase is energy 
independent. This may happen both in the singlet as well 
as in the triplet channels with j — I. The short distance 
phase is then fixed by reproducing the scattering length 
in the singlet channel and one of the three scattering 
lengths in the triplet channel. In the case where one has 
two negative, i.e. attractive, eigenvalues (this can only 
happen in triplet channels) there are two undetermined 
short distance phases which can be fixed by using the 
corresponding three scattering lengths. The case of two 
positive, i.e. repulsive, eigenvalues does not allow to fix 
any scattering length. The case with two different signs 
for the eigenvalues fixes one scattering length only. Note 
that in this construction and for two coupled channels 
there is no intermediate situation where the solution is 
specified by just two scattering lengths; one has either 
zero, one or three. 



Although our arguments are entirely based on analyti- 
cal calculations, one should mention that our conclusions 
are in agreement with the findings of Ref. ||42j for the 
OPE case. There, counterterms beyond the ones dictated 
by Weinberg's power counting are included in the ^Pq, 
— and '^D2 waves to ensure renormalizability on 
numerical grounds. As we will see below, our renormal- 
ized phase shifts for special case OPE reproduce essen- 
tially their results, although our TPE non-perturbatively 
renormalized amplitudes go beyond these results. 



Another issue is that of the establishment of a theoret- 
ically compelling and mathematically consistent power 
counting which also provides phcnomenological success. 
This has been the goal of much the EFT activity in recent 
years. Despite the fact that our OPE is mathematically 
identical to the one in Ref. jl^l where a strong empha- 
sis on power counting has been made, our motivation 
is slightly different. Actually, these authors argue that 
a consistent scheme for TPE might be achieved within 
a perturbative framework, using the non-perturbative 
OPE distorted amplitudes as the leading order approx- 
imation. This is theoretically appealing and the issue 
was thoroughly discussed within the coordinate space ap- 
proach in our previous paper on the central waves |43| . 
There, it was pointed out that with enough countert- 
erms such a program could be pursued although orthog- 
onality was violated and results did not exhibit a clear 
improvement as compared to the fully iterated poten- 
tials. The reason was the appearance of non-analytical 
dependences on the would-be dimensional power count- 
ing parameter, a situation that has not been foreseen in 
the standard EFT set up. This suggests that discussion 
on power counting and the systematics of EFT is not 
yet over. Therefore, and as we did in our previous work, 
we focus more on establishing long range model indepen- 
dent correlations, leaving the possible establishment of a 
satisfactory power counting for future studies. 
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C. Regularization methods 

In principle, it is possible to implement the short dis- 
tance behaviour of the wave functions, Eq. (|23Jl, if one 
goes to sufficiently small distances, or if the short dis- 
tance behaviour of the the wave function is improved |4lj | . 
Computationally, the implementation of short distance 
regulators is mostly straightforward. The attractive or 
repulsive nature of the potentials at short distances re- 
quires different choices of regulators [4ll |43| . For a one- 
channel repulsive singular potential we use the regulator 

<M^[±1, (27) 
Uk[a) a 

This condition ensures orthogonality of wave functions 
with different energy. For the attractive singular case, 
we integrate in from infinity at zero energy down to a 
given boundary radius, a, impose orthogonality at the 
boundary by matching logarithmic derivatives 

<(«) /9<5^ 
Uk{a) uo[a) 

and then integrate out at finite energy. In the coupled 
channel case we extend the method by applying the one 
channel regularization to the short distance eigen func- 
tions, Eq. 

D. Fixing of parameters and renormalization 
conditions 

Fixing of the short distance phases requires some 
renormalization condition. As we have said, an appealing 
choice is to impose this condition at zero energy. The way 
to proceed in practice is quite straightforward, although 
tedious given the large number (27) of partial waves con- 
sidered in this work. In the singlet channel case and 
for an attractive short distance singularity, one starts at 
zero energy and integrates in from large distances ~ 15fm 
with a given scattering length until a short boundary ra- 
dius ~ O.lfm. At finite energy one integrates out match- 
ing the wave function to the zero energy solution at the 
short distance boundary generating a phase shift from 
a given prescribed scattering length. Of course, in this 
method one has to check for cut-off independence (taking 
r = 0.1 — 0.2fm proves enough). For the coupled chan- 
nel case one proceeds along similar lines and the proce- 
dure has been described in great detail in our previous 
works pH E^ l for the j = 1 channel. The method relies 
heavily on the superposition principle of boundary con- 
ditions and we use here the extension of that method to 
higher partial waves. One of the advantages of our ap- 
proach is that we rarely have to make a fit to the data; 
any phase shift has by construction the right threshold 
behaviour in the case where the potential at short dis- 
tances is attractive. For the repulsive potential case the 
scattering length is predicted entirely from the potential. 



TABLE IL The number of independent parameters for differ- 
ent orders of approximation of the potential. The scattering 
lengths are in fm'+' +1 and are taken from Nijmll and Reid93 
potentials ^3 in Ref. 13 . We use the (SYM-nucIear bar) 
convention, Eq. 1161 . 



Wave 


a Nijmll (Reid93) 


LO 


NLO 


NNLO 


'So 


-23.727(-23.735) 


Input 


Input 


Input 




-2.468(-2.469) 


Input 


— 


Input 


'Pi 


2.797(2.736) 


— 


— 


— 


''Pi 


1.529(1.530) 


— 


Input 


Input 


'Si 


5.418(5.422) 


Input 





Input 




6.505(6.453) 








Input 


El 


1.647(1.645) 


— 


— 


Input 




-1.389(-1.377) 


— 


Input 


Input 




-7.405(-7.411) 


Input 


Input 


Input 


'P2 


-0.2844(-0.2892) 


Input 


Input 





'F2 


-0.9763(-0.9698) 











E2 


1.609(1.600) 


— 


— 


— 


'F, 


8.383(8.365) 


— 


— 


— 


'F, 


2.703(2.686) 


— 


Input 


Input 


'D, 


-0.1449(-0.1770) 


Input 


— 


Input 




4.880(4.874) 


— 


— 


Input 


Es 


-9.695(-9.683) 


— 


— 


Input 


'G, 


-3.229(-3.210) 




Input 


Input 


'G, 


-19.17(-19.14) 


Input 


Input 


Input 


'F, 


-0.01045(-0.01053) 


Input 


Input 




'H, 


-1.250(-1.240) 








Ei 


3.609(3.586) 








'H, 


28.61(28.57) 








'Hs 


6.128(6.082) 




Input 


Input 


'G, 


-0.0090(-0.010) 


Input 




Input 


'h 


10.68(10.66) 






Input 


E5 


-31.34(-31.29) 






Input 



In any case, discrepancies with the data can be attributed 
to the potential. 

Inspection of Table |n] illustrates the situation for the 
LO, NLO, and NNLO approximations to the potential. 
We show the scattering lengths in all partial waves as 
determined in our previous work j48l | together with the 
corresponding eigenvalues for the leading short distance 
coefficients in the LO (OPE), NLO and NNLO approx- 
imations to the potential. In the NNLO one must also 
specify the values of the chiral constants ci, C3 and C4. 
We use for definiteness the values of Ref. ^ , since as we 
saw in Ref. ,43] they provide a reasonable description of 
deuteron properties. 

E. Details on the numerical procedure 

The integration of the coupled differential equations re- 
quires some care, particularly in the vicinity of the short 
distance singularities. In the case of attractive singular- 
ities due to the increasing oscillations the wave function 
has to be sampled with great detail at a rate similar to 
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the size of the oscillations. For the repulsive case, one 
must stop at sufficiently large distances due to the ex- 
ponential suppression of the wave function. Another im- 
portant condition has to do with preservation of in and 
out reversibility of the integration. This last require- 
ment guarantees that for attractive channels, where the 
scattering length is supplied as an input parameter, the 
threshold behaviour of the phase shift is consistent with 
that given scattering length. 

Another problem one has to face for high partial waves 
is related to the practical influence of the scattering 
length on the calculated phase shifts. In principle, and 
for an attractive singular potential, the scattering length 
needs to be specified. For the one channel case, this is 
done by integrating in the zero energy large distance so- 
lution, valid for r ^ 2/m.„ 

u{r) r-' , (29) 

The long distance irregular solution dominates, unless 
ai is anomalously large , i.e. a/(TO7r/2)''^'+^^ 3> 1 , so 
that when integrating in much of the regular solution 
will be lost and the result will be rather insensitive to 
the value of ai provided it is of normal size. This fact 
becomes relevant in the numerical calculations if the long 
distance cut-off is taken to be exceedingly large. To avoid 
this situation we take typically i?max — 15fm for large I. 



III. RESULTS FOR THE PHASE SHIFTS 

A. Numerical parameters 

For our numerical calculations we take /tt — 92.4MeV, 
TO = 138.03MeV, 2^„p = M = MpMr^/{Mp + M„) = 
938.918MeV, gA = 1.29 in the OPE piece to account 
for the Goldberger-Treimann discrepancy and gA — 1.26 
in the TPE piece of the potential. The correspond- 
ing pion nucleon coupling constant takes then the value 
girNN = 13.083 (i.e. gA = 1.29) according to the Ni- 
jmegen phase shift analysis of NN scattering [s^l- The 
values of the coefficients ci, C3 and C4 used along this 
paper can be looked up in Tabled for completeness. The 
potentials in configuration space used in this paper are 
exactly those provided in Ref. [1,^1^ but disregarding 
relativistic corrections, M/E — > 1 The potentials are 
listed in Appendix ^ for completeness. The short dis- 
tance van der Waals coefficients for all channels studied 
in the present work are presented in Appendix [CI The 
output of such a channel by channel analysis is briefly 



As mentioned in our previous work l4,?l these effects are tiny 
for the deuteron. For central waves they are about 0.2 ^ at the 
maximum CM momentum p = 400MeV. This trend is general 
also for peripheral waves. 



summarized in Tabled where we indicate which scatter- 
ing lengths are used as input parameters according to the 
discussion given in Sect. Low energy parameters for 
the high quality potentials ^t'I have been obtained 
in Ref. 48] . We will use the Nijm II values, but to give 
an idea on the expected lower uncertainties on those pa- 
rameters we also list the Reid93 values. Probably the 
real uncertainties are much larger since the actual value 
of these low energy parameter will depend upon which 
long range physics is included in the high quality poten- 
tials where explicit TPE effects have not been included, 
as we do in the present work ^. 



B. The deuteron channel revisited 



Before embarking in the full fledged discussion of all 
partial waves, it is interesting to reanalyze first the 
'^Si Di channel already studied in our previous work 
on the deuteron 0,01 ■ There, we used the orthogonal- 
ity to the deuteron bound state. The scattering lengths 
cko2 — 1.67 fm^ and a2 — 6.6 fm^ were deduced from the 
experimental deuteron binding energy, the asymptotic 
D/S ratio and the S-wave scattering length ao. These 
values turned out to be a bit off the values deduced from 
the Nijmll and Reid93 potentials ^ (see Table|lll). Nev- 
ertheless, the intermediate energy region turned out to 
be better described than the low energy behaviour sug- 
gested. In the present work we choose instead to build 
scattering states which are orthogonal to the zero energy 
states, so deuteron properties can be deduced, as done in 
table im In Fig. ^ we show the results when either the 
zero energy or the deuteron bound state are used as ref- 
erence states. One obvious lesson from this comparison 
is that phase shifts, particularly the Ei channel, may be 
better described in the intermediate energy region if the 
deuteron is used as a reference state, despite the fact that 
the threshold behaviour is a bit off. This is maybe ex- 
plained by the observation that ao2 and 02 encode higher 
energy information about the system than ao or 7 ^, so 
the latter parameters are more suited to obtain an effec- 
tive description of the system. This feature will become 
evident in other partial waves. 



^ This will generate slight inconsistencies in the TPE results of 
Sect. 1111 D 31 which will be amended by a small modification of 
the threshold parameters, yet larger than the discrepancies be- 
tween the threshold parameters for Nijmll and Reid93 potentials 
obtained in Ref. l48l . 

^ It should be note that ao2 and a2 are related with the behaviour 
of the scattering amplitude at order and fc* respectively, rel- 
ative to 00 . 
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TABLE III: Deuteron properties for the OPE and TPE potentials. OPE(O) refers to the deuteron computed by orthogonality 
to the zero energy scattering states, fixing ao to its experimental value, while in OPE(B) the computation is made by fixing 
7 to its experimental value, constructing the corresponding bound state. Similarly, TPE(O) refers to the deuteron computed 
by orthogonality to zero energy states, fixing ao to its experimental value and ao2 and 02 to their Nijmegen II values, while in 
TPE(B) the computation is made by fixing 7, -q and ao to their experimental values. The errors quoted in OPE(O) corresponds 
to the uncertainty in the value of the scattering length, while in OPE(B) the errors correspond to changing the cutoff in the 
0.1 — 0.2 fm range. The errors quoted in both TPE computations reflect the uncertainty in the non-potential parameters 7, rj 
and ao only. We take set IV for the LEC's in the TPE calculation. Experimental values can be traced from Ref. 57]. 



Set 


7 (fm-^) 


V 


As (fm-i/^) 


Tm (fm) 


Qd(fni^) 


Pd (%) 


ao (fm) 


ao2(fm^) 


a2(fm^) 


OPE(O) 


0.2274(4) 


0.02564(4) 


0.8568(10) 


1.964(3) 


0.2796(3) 


7.208(12) 


Input 


1.754(7) 


6.770(7) 


OPE(B) 


Input 


0.02633 


0.8681(1) 


1.9351(5) 


0.2762(1) 


7.31(1) 


5.335(1) 


1.673(1) 


6.693(1) 


TPE(O) 
TPE(B) 


0.2322(3) 
Input 


0.02531(9) 
Input 


0.8891(4) 
0.884(4) 


1.968(3) 
1.967(6) 


0.2723(3) 
0.276(3) 


7.24(13) 
8(1) 


Input 
Input 


Nijmll 
1.67(4) 


Nijmll 
6.6(4) 


Nijmll 
Reid93 


0.231605 
0.231605 


0.02521 
0.02514 


0.8845 
0.8845 


1.9675 
1.9686 


0.2707 
0.2703 


5.635 
5.699 


5.418 
5.422 


1.647 
1.645 


6.505 
6.453 


Exp. 


0.231605 


0.0256(4) 


0.8846(9) 


1.971(6) 


0.2859(3) 




5.419(7) 








ElabIMbVI ElabIMsV] ELABlMeV] 



FIG. 1: Dependence of the ^Si — ^ Di channel (SYM-nuclear bar) phase shifts for the NNLO TPE potential on the reference 
state used to orthogonalize the scattering state compared to the corresponding phases of the database of Ref. |4^ li^. Label 
TPE(O) means the zero energy reference state with ao = 5.418fm, ao2 = 1.647fm'^ and a2 — 6.505fm^. Label TPE(B)stands for 
the deuteron bound reference state with the experimental binding energy, asymptotic D/S ratio together with ao = 5.418fm 
(corresponding to ao2 = 1.67 fm^ and a2 = 6.6 fm^ 01-) 



C. Cut-ofF dependence 

In Figs. 121 El E El El and [7\ we show the results of our 
calculation for all partial waves with j < 5 as a function 
of the nucleon LAB energy. For definiteness we use the 
chiral constants ci , C3 and C4 of Ref. P| (Set IV) which 
already provided a good description of deuteron prop- 
erties after renormalization at NNLO. This choice 
allows a more straightforward comparison to the N'^LO 
calculation of Ref. ^] with finite cut-offs. Unless other- 
wise stated, the needed low energy parameters for these 
figures are always taken to be those of Ref. for the 
Nijmll potential (see Table Hljl . 

In order to test the stability of the phase-shifts against 
changes in the short distance cut-off parameter, we 
show in Figs. El El El El El and [3 and similarly to the 
OPE study in momentum space of Ref. > the cut-off 
dependence for fixed values of the lab energy both for 
the OPE as well as for the TPE potentials. This is done 
in the range 0.15fm < Rs < 1.5fm. If we identify this 
short distance cut-off with the sharp momentum cut-off 



A = tt/2Rs 45], the smallest boundary radius, ^ 0.15fm, 
corresponds to a maximum cut-off A ^ 2GeV. This is 
much larger than the cut-offs used in Refs. H, S S IS 
0, ITo| but comparable to the exponential cut-off used in 
Ref. Hy\ for the renormalization of the OPE potential ^. 
Note that the limit Rs may be taken independently 
for any different channel. 

The evolution of the increasingly oscillating wave func- 
tion in the attractive case can be identified with the 
cycles (improperly called limit-cycles, see footnote 5 in 



^ There, a cut-off has been introduced according to the rule in the 
potential V{k',k) — » e"*" V{k'k)e~'" and counterterms 
have been added. To get an order of magnitude of the equivalent 
sharp cut-off A we estimate the linear divergence at zero energy 
in the contact theory, 

K = r = = 0.762A , 

and also using A = n/2Rs we get A = l/(0.48i?s) . 
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FIG. 2: np (SYM-nuclear bar) Phase shifts for the total angular momentum j — 0. OPE (left) and chiral TPE (middle) as 
a function of the cut-off radius Rs for fixed LAB energies, Slab ~ 10, 50, lOOMeV. OPE and chiral TPE (left) renormalized 
(i.e. Rs — > 0) phase shifts as a function of the LAB energy compared to the Nijmegen partial wave analysis |46L li^ . 



Ref. liol) described in Refs. 'sl S E by looking 
at suitable logarithmic combinations of the wave func- 
tions. The cycles documented in Ref. ^3 momentum 
space can be mapped into the coordinate space cycles by 
relating the coordinate and momentum space cut-offs. 

Generally speaking, the inclusion of chiral TPE effects 
generates smoother limits as compared to the OPE re- 
sults, as one would expect. We have checked that for 
short distance repulsive (eigen)channels results are not 
very sensitive to the choice of the regulator for small val- 
ues of Rs- As we also see from the figures, the conver- 
gence depends both on the partial wave as well as on the 
energy. As expected, the needed value of the short dis- 
tance cut-off Rs for which stability is achieved is rather 
high for peripheral waves, Rs ~ Another feature 

of the calculation are the observed stability plateaus for 
a number of partial waves. This trend has also been 
noted in previous works with finite cut-offs [Tlj where 
there appear sequential cut-off windows. In coordinate 
space this is originated by the almost self-similar pattern 
of the short distance oscillations of the wave function 
which suggest a sequential and faster convergence mod- 
ulo cycles 4Qj. 

Let us remark at this point that the existence of an 
Rs — *■ limit does not necessarily mean a plateau-like 
approach to it. This is the case, for example, of the ^Sq 
wave, which for OPE shows a linear dependence on the 
cut-off due to the mild 1/r singularity of the potential, 
generating a linear-like behaviour which corresponds to 
the ratio of regular (~ r) and irregular (~ 1) solutions 



at the origin ®. A similar behaviour can be found on 
other singlet waves in which the OPE potential also be- 
haves as 1/r, but highly attenuated by the influence of 
the centrifugal barrier. 

Finally, let us note that there are some channels where 
the phase shifts exhibit a very strong dependence on the 
regulator ^. 

D. Renormalized phase shifts 

1. LO (OPE) 

In Figs. 13 Em El El and □ we also compare the OPE 
(LO), the NNLO TPE and the Nijmegen phase shift anal- 
ysis |4a,|43. As noted in Tab. Ullin some cases with at- 
tractive singular potentials some scattering lengths must 



Anyway, the lack of a clear plateau in this wave becomes obvious 
in the coordinate space treatment. Assuming the relationship 
A ~ l/(0.48-Rs) for a gaussian cut-off ( see footnote 171 between 
the momentum and coordinate space cut-offs, a linear depen- 
dence of the phase shifts on the Rg coordinate cut-off would map 
into a 1/A dependence in momentum space, which might be re- 
garded as a plateau in a sufficiently thin cut-off window.. Note 
that going from Rb = 0.2fm to Rs = O.lfm corresponds equiv- 
alently to double the momentum space cut-off from ~ 2GeV to 
~ 4GeV. 

The jump in the evolution of the OPE potential in the chan- 
nel around Rg = 0.3fm, Fig. ISlresembles a coupled channel reso- 
nance, corresponding to tunneling across the centrifugal barrier 
into the short distance attractive singularity. 
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FIG. 3: Same as Fig.Hbut for j = 
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FIG. 4: Same as Fig.Hbut for j = 2. 
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FIG. 5: Same as Fig.Hbut for j = 



3. 
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FIG. 6; Same as Fig.Hbut for j = 



4. 




FIG. 7: Same as Fig.Hbut for j = 5. 
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be specified, in order to determine the phase shifts, but 
for repulsive singular potentials the scattering lengths 
and hence the phase shifts are fully determined from the 
potential. In the coupled channel case where only one pa- 
rameter should be fixed we have chosen, as indicated in 
Table Hll to take the scattering length of the correspond- 
ing partial wave with the lower orbital angular momen- 
tum. As we see from Figs. El El H El El and OPE does 
a relatively good job for the phases when compared to 
the Nijmll results, up to a reasonable energy. This cal- 
culation extends our previous results ^| using the same 
regularization for the singlet ^Sq and triplet '^Si Di 
channels. 

The LO results corresponding to static OPE potential 
have also been obtained recently in momentum space by a 
solution of the Lippmann-Schwinger equation in Ref. |42l | 
for j < 3. These authors see that in the limit A ^ oo 
(in practice A = 4GeV) it is always possible to adjust a 
counterterm in such a way that the phase shifts are cut- 
off independent. They also find that the needed coun- 
terterm does not correspond to the expectations based 
on Weinberg's dimensional power counting argument, so 
that one is forced to promote counterterms which are 
of higher order in Weinberg's counting to make the the- 
ory free of short distance ambiguities. This proposal not 
only fits quite naturally into our analysis of short dis- 
tance boundary conditions, but can also be anticipated 
by just looking at the short distance behaviour of the 
potential. In general, we reproduce their results for the 
phase shifts using our boundary condition regularization 
(our shortest distance cut-off is typically a — O.lfm for 
OPE). This is precisely one of the points of renormaliza- 
tion; different regularization methods should yield iden- 
tical results when the regulator is removed provided the 
same renormalization conditions are imposed. Note that 
in our case whenever a scattering length must be provided 
we exactly construct the phase shift as to reproduce the 
threshold behaviour of the Nijmegen phases [461 It^ by 
exactly fixing the scattering length (the renormalization 
condition). This requires solving the zero energy prob- 
lem by integrating in with the given scattering length, 
and matching at short distances the finite energy prob- 
lem to finally determine the phase shift by integrating 
out. In this approach we never make a fit. In the ap- 
proach of Ref. |43] counterterms are adjusted to fit the 
phases in the region around threshold. Although this is 
in spirit the same renormalization condition to fix the 
counterterms, we expect some numerical discrepancies, 
due to the fact that the threshold parameters in Ref. 
may be slightly different to ours. 



distance behaviour of the NLO potential corresponds to 
l/r^ repulsive eigenpotentials. This feature explains the 
relatively small maximal cut-offs allowed in NLO calcu- 
lations in momentum space. As stressed on our previous 
work there are at least two scenarios where the prob- 
lem may be overcome. One possibility appeals to the 
role of the A resonance and the fact that its contribution 
to C3 and C4 scales as the inverse of the iVA splitting 
A ~ 2m7r as found in Ref. 0, flsL [56| . In the A count- 
ing the C3 and C4 contributions to the NNLO deltaless 
potential become actually NLO contributions, and the 
short distance behaviour becomes a 1/r^ attractive sin- 
gularity. The second scenario has to do with the influence 
of relativity beyond a truncated heavy baryon expansion, 
since according to Refs. [23| one has a relativis- 

tic l/r^ van der Waals short distance behaviour with 
attractive-repulsive eigen potential meaning that as in 
the OPE case one has one free parameter. Calculations 
taking into account these effects in all partial waves are 
currently underway |52|. 

3. NNLO (chiral-TPE) 

We turn now to the NNLO calculation which are the 
genuine predictions of ChPT because they contain the 
chiral constants ci, C3 and C4 (see e.g. Table and for 
definiteness we will use mainly Set IV Q in our analy- 
sis . Results for the TPE renormalized phase shifts are 
presented in Figs. El El E El El a-nd [7| Some expected 
features do indeed occur. Peripheral waves are slightly 
modified by going from OPE to the chiral TPE potential. 
On the other hand, low partial waves are also improved 
in the low energy region. For instance, the ^Sq phase 
has an attractive singular interaction, requiring fixing the 
scattering length. The difference in the curves is mainly 
related to the difference in the effective r ang e which im- 
proves when going from OPE to TPE [43 ■ This is a 
rather general feature, the error at low energies is con- 
trolled by the low energy threshold parameters, like the 
effective range and others. If one looks at the '^Pq chan- 
nel, we see that there is improvement but not as dramatic 
as in the ^Sq channel. 

As we have said, in singular repulsive channels, which 
at NNLO correspond to the ^Pi, ^Pa and ^Hr^ singlet 
states, and to the — and '^Pj — ^H^ triplet states, 
the phase shift and also the scattering length are entirely 
determined by the potential. So, these phases are a good 
place to study the influence of different values for the chi- 
ral constants, ci, C3 and C4, presented in Tabled In Fig. El 
we show this dependence for these special partial waves. 



2. NLO (TPE) 



Regarding NLO we do not show the results as they 
fail completely to describe the data in triplet ^Si Di 
channel. The problem we found already "431 in the triplet 
^Si Di channel persists in other channels; the short 



By genuine we mean that the NNLO potential contains parame- 
ters which are relating irN and A'^A'^ data in some intricate way. 
The fact that we are using the parameter Set IV Q is because it 
nicely reproduces the deuteron properties. One could, of course 
improve on this by a large scale fit to the data 
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FIG. 8: Dependence of some phases for the chiral TPE potential on the scattering lengths compared to the Nijmll phases l4fiLl47l|. 
We use Set IV of chiral constants. Label TPE (anijm) means scattering lengths of TablelTTlare taken as deduced in Ref. pal , while 
TPE (Qtuncd) stands for the values tuned to fit the phases in the intermediate energy region. For ^Po we take a\ = — 2.670fm^, 
for ^Pi we take ai = 1.692fm^ and for ^D2 we take 0:2 = -1.666fm^ 



As we see, the ^Pi phase exhibits a strong dependence on 
the parameter set, while -^Fs and ^H^ are less sensitive 
to this particular choice. The strong dependence in the 
^Pi channel suggests that this may be an ideal place to 
fit the chiral constants, since the scattering lengths are 
fixed. We will not attempt such determination of the 
chiral constants here because that would require realistic 
error estimates of the phase shifts. 

If we restrict to the spin singlet channels we see that 
there is very good agreement for higher peripheral waves, 
^i?5, ^Gi and ^^3. This is expected from perturbative 
calculations. Note, however, that unlike perturbation 
theory we fix by construction the scattering lengths for 
the case of singular and attractive potentials. Some inter- 
mediate waves, such as ^D2, which potential is singular 
attractive, are badly reproduced despite the fact that the 
threshold behaviour is in theory reproduced since we use 
the corresponding scattering length as input. Actually, 
for these waves the TPE result seems to worsen the OPE 
prediction. Presumably this is an indication either on 
the inadequacy of the (Nijmll) scattering lengths used 
as input for NNLO or on the importance of N'^LO con- 
tributions. Let us note that the Nijmll potential does 
not incorporate explicit TPE effects in their long range 
part. In fact, if we take a slightly different scattering 
length, 02 = — 1.666fm^ instead of the values deduced 
in Ref. Hg, = -1.389fm^ for the Nijmll potential 
a rather good agreement with the Nijmegcn analysis is 
obtained for the phase shift (see similar results for 
^Po and ^Pi waves in Fig. ISJ. Although the small dif- 
ference between the fitted and experimental values for 
the scattering length could also be explained by N^LO 
corrections, suggesting that they are not large, a definite 
conclusion cannot be drawn in the absence of a large scale 
fit 



These general trends are confirmed in the triplet chan- 
nels, where in high partial waves there is an overall im- 
provement when going from OPE to TPE. In some cases, 
like in the , £2, and ■^P^ the improvement is 
rather satisfactory all over the energy range. However, 
the theory has notorious problems in the "^Pi and ei and 
to a lesser extent in the ^D^ and Ei channels if one insists 
on keeping the scattering lengths of the Nijmll potential. 
As before, small changes in the scattering lengths allow 
for an overall improved description as can be deduced 
from Fig. |S1 in some particular cases (see also Fig^. 
This suggests that higher orders in the potential may be 
needed. This fact was pointed out in our previous work 
on the central phases, where the NNLO potential made 
almost the effective range although there was a statisti- 
cally significant discrepancy to the experimental number, 
which called for the inclusion of N^LO terms. This may 
possibly happen also in some higher partial waves and it 
would be interesting to see whether improved long dis- 
tance potentials might account for the observed discrep- 
ancies to phase shifts provided the scattering lengths are 
kept to their physical values. 

As we have shown (see Fig. ISjl, small changes in the 
scattering lengths indeed allow for a better description 
of the phases in the intermediate energy region. On the 
other hand, we would expect our description to become 
increasingly better for lower energies. This situation is 
a bit disconcerting. Given the similarity between the 
scattering lengths computed in Ref. 48] for the Nijmll 
and Reid93 potentials it seems unlikely that potential 
models yield a completely off value for the a's in non- 
central waves, but one must admit that errors will in 
general be larger than the difference between these two 
potential values suggests as already argued above. If one 
takes into account the fact that both potentials include 
similar long range physics, this means that the true error 



This also applies to the non-static OPE corrections which ac- 
count for about 0.1° at Elab = 200MeV. The effect can be 
mocked up by even tinier readjustments of both the scattering 



lengths as well as the chiral couplings ci , C3 and C4 than deduced 
from inaccuracies in the Nijmll potentials. 
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FIG. 9: Dependence of the NNLO (SYM-nuclear bar) phases on the chiral constants. These are the only channels where the 
potential is singular repulsive. 



could be larger due to systematic uncertainties in their 
short range part (non OPE). 

Nevertheless, let us mention that current calculations 
involving chiral potentials not only ignore this possible 
disagreement at threshold but they in fact modify the 
corresponding scattering lengths since the counterterms 
are determined by a fit to the phase shifts in the re- 
gion above threshold with no obvious control on the low 
energy parameters (see e.g. Ref. "42]). The arguments 
above do not prove that taking slightly different scatter- 
ing lengths than those suggested by the high quality po- 
tentials is a legitimate operation but at least shows that 
no more assumptions are made. From this viewpoint it 
might be profitable to study the impact on those calcu- 
lations of either imposing exact threshold behaviour or 
alternatively evaluating the threshold parameters them- 
selves. 



E. Remarks on the perturbative nature of 
peripheral waves 

The numerical coincidence of our non-perturbative cal- 
culations with perturbation theory expectations 0, |3 , 
although quite natural on physical grounds, deserves 
some explanation on the basis of the formalism and the 
relevance of short distance singularities. Indeed, the at- 
tractive character of the singular NNLO potentials at the 
origin implies a non-trivial boundary condition of the 
form of Eq. (|24|1 . which cannot be reproduced to any 
given order in perturbation theory, at least without the 
inclusion of extra counterterms in the perturbative ex- 
pansion, a point which will be further discussed at the 
end of this section. This point was previously illustrated 
in Ref. .36] for s- waves and also in our previous work on 
the renormalization of the OPE 41] by comparing the ex- 
act deuteron wave functions with the perturbative ones. 
There, one observes that the first order perturbative cal- 
culation provides finite results, but the expansion at sec- 
ond order produces divergent results due to the short 
distance non-normalizable D— wave component. Thus, 
observables cannot, strictly speaking, be analytical func- 
tions of the coupling (for the purpose of discussion we 



could visualize the problem by thinking of singularities 
of the sort + g'^logg'^). This does not mean that for 
the physical range of couplings the non-analytical contri- 
bution is necessarily large numerically. For instance, in 
the deuteron channel the residual non-analytical higher 
order terms happens to be numerically sizeable even for 
a weakly bound deuteron. 

Based on the results of Ref. there is no reason 

to expect that higher partial waves will not exhibit this 
failure of perturbation theory at some finite order. Nev- 
ertheless, the perturbative short distance behaviour of 
higher partial waves tames the singularity due to the 
kinematical suppression. This is a perturbative long 
distance feature where the centrifugal barrier dominates. 
The point is that this short distance behaviour is not in- 
variant order by order in strict perturbation theory for a 
singular potential and, actually, one finds a short distance 
enhancement of the wave function even in perturbation 
theory. So, one expects that the perturbation theory on 
a singular potential will diverge at some finite order also 
for high partial waves. In Appendix^we show that this 
is indeed the case; for a singular potential diverging like 
l/r" (n > 2) and a partial wave with angular momen- 
tum the perturbative expansion diverges at A:— th order 
in perturbation theory provided fc > (2^ + l)/(?i — 2). 
This estimate provides the order at which, if desired, 
a long distance perturbation theory on boundary con- 
ditions might be applied as discussed previously for the 
deuteron channel [4l| . Using the techniques developed in 
Ref. ^3 to make perturbation theory on distorted OPE 
central waves it would be interesting to see, as claimed 
by renormalization arguments on the OPE |54l |. whether 
such an expansion is indeed possible. 

Having established that perturbation theory will di- 
verge at some finite order, we would like now to under- 
stand why it still can accurately represent the full non- 
perturbative solutions obtained numerically. The reason 
can be found in the very efficient way how the short dis- 
tance singularity of the potential makes short distances 
to be inessential in the wave function for the regular non 
perturbative solution. For high angular momenta and at- 
tractive singular potential the wave function senses the 
singularity after tunneling through the barrier, an expo- 
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nentially suppressed effect. In perturbation theory the 
effect is just substituted by the core provided by the cen- 
trifugal barrier. 



IV. CONCLUSIONS 

In the present paper we have analyzed the renormal- 
ization of non-central waves for NN scattering for the 
OPE and chiral TPE potentials. This calculation extends 
our previous studies on central phases and the deuteron 
for OPE and TPE potentials presented in Refs. [U El] 
respectively. As already stressed in those works the re- 
quirement of finiteness of the scattering amplitude as well 
as the orthogonality of wave functions impose tight con- 
straints on the allowed structure of counterterms for a 
given potential. Using the standard Weinberg count- 
ing for the potential, the counterterm structure is de- 
duced and does not generally coincide with the naive 
expectations. In some cases forbidden counterterms in 
the Weinberg counting must be allowed whereas 
in some other cases allowed counterterms must be ex- 
cluded ^3 ■ Finite cut-off calculations based on the Wein- 
berg counting allow to introduce counterterms which are 
usually readjusted to globally fit the data but are for- 
bidden by finiteness and orthogonality, in renormalized 
calculations. The success of the original counting relies 
heavily on keeping finite the cut-off, while at the same 
time it is usually emphasized that low energy physics 
does not depend crucially on short distance details. As 
we have argued, these two facts are mutually contra- 
dicting; the standard Weinberg counting is incompati- 
ble with exact renormalization, i.e. removing the cut-off, 
as was suggested in Re f. [ 55l within a perturbative set 
up and shown in Ref. '4^ non-perturbatively, at least 
in the heavy baryon expansion and when only nucleons 
and pions are taken into account. This feature changes 
when relativistic effects and A degrees of freedom are 
taken into account, showing that perhaps renormaliza- 
tion, i.e. the independence on short distance details may 
be a strong condition on admissible potentials. In this 
regard we find that, as one would expect, the cut-off de- 
pendence is milder for the chiral TPE potential than for 
OPE potential. This suggests that higher order correc- 
tions become even more cut-off independent. Indeed, the 
finite cut-off N^^LO calculations of Ref. jTlj do exhibit 
this feature in spite of the strong cut-off dependence ob- 
served at lower orders. 

Using this modified Weinberg counting, the quality of 
the agreement and improvement depends on the partic- 
ular partial wave. High partial peripheral waves, when 
treated non-perturbatively reproduce the data fairly well 
and deviations from OPE to TPE are small, as one would 
expect in a perturbative treatment. Nevertheless, we 
have also shown that regardless on the orbital angular 
momentum, there is always a limit to the order in pertur- 
bation theory for which finite results are obtained. The 
divergence is related to an indiscriminate use of the per- 



turbative expansion, and not to an intrinsic deficiency in 
the definition of the scattering amplitude. Thus, also for 
peripheral waves the phase-shifts are perturbatively non 
renormalizable while they are non-perturbatively renor- 
malizable. This result extends a similar observation for 
the deuteron [40ll4lj|. Nevertheless, we have also argued 
why convergent perturbative calculations to finite order 
are useful and may even provide accurate descriptions 
when compared to the non-perturbative result. 

Unlike naive expectations, it is not always true that 
after renormalization the NNLO TPE phases improve 
over OPE ones if one insists on keeping the scattering 
lengths required by finiteness to the same physical values 
as those extracted ij48j from the high quality Nijmegen 
potentials|47j. This renormalization condition at zero en- 
ergy has been adopted to highlight the difference between 
these potentials j43 and the chiral NNLO singular poten- 
tials jJI . Remarkably, using zero energy to fix the param- 
eters has never been considered before within the chiral 
potentials approach to NN scattering, thus some of the 
problems we find and discuss have not even been identi- 
fied so far. Actually, we find that some partial waves such 
as ^I?2 and ^Pi are particularly sensitive to the value of 
the scattering length. In fact, it is found that small de- 
viations of the scattering lengths at the few percent level 
in these partial waves improve dramatically the descrip- 
tion in the intermediate energy region. The improvement 
can also be achieved in other partial waves by suitably 
tuning the scattering lengths in all the channels charac- 
terized by singular attractive interactions. This means 
that the absolute error is small up to i?LAB ^ lOOMeV. 
Three pion exchange effects should become relevant at 
about CM momentum of fc = 3m/2 which corresponds 
approximately to this LAB energy. The modification cor- 
responds to change the renormalization condition to some 
finite energy, or maximizing the overlapp between the 
chiral phase shifts and the fitted ones in a given energy 
window, very much along the lines pursued in previous 
works. However, changing the scattering lengths pro- 
duces large relative errors near the threshold. At this 
point the discussion on errors on the phase shifts be- 
comes a crucial matter, particularly in the low energy re- 
gion. In this regard, it seems likely that the difference in 
low energy threshold parameters determined in Ref. [4^ 
for the Reid93 and Nijmll in all partial waves with j < 5 
provides a lower bound for the true error. Obviously, a 
meticulous error analysis of these threshold parameters 
would be very helpful. 

We have also found that some partial waves, with re- 
pulsive singular interactions and where no free scatter- 
ing lengths are allowed, are particularly sensitive to the 
choice of chiral constants ci, C3 and C4. This suggests 
that a fit of the chiral constants to these partial waves 
may be possible. To do so, and again, a realistic estimate 
on the errors of the phase shifts would be mandatory. 
According to our findings on the deuteron for the chiral 
TPE potential it is quite likely that, if such error es- 
timate was reliably done, theoretical determinations for 
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deuteron observables with unprecedented precision based 
on chiral potentials might be achieved. This issue is be- 
ing curre ntly under consideration and is left for future 
research 53]. 

From a practical viewpoint there is a potential disad- 
vantage in requiring exact renormalization for the ap- 
proximated long distance chiral potentials, due to the 
tight constraints imposed by finiteness on the short dis- 
tance behaviour of the wave functions. To some extent, 
although the chiral potentials are motivated by the Ef- 
fective Field Theory idea, these additional conditions re- 
mind also aspects of renormalization of fundamental the- 
ories. This is not entirely surprising since we expect the 
chirally based potentials to resemble the true NN po- 
tential, at least at sufficiently long distances. For in- 
stance OPE is a true long distance contribution. Full 
TPE would also be a true long distance part, which is 
known in an approximate manner within the current 
ChPT schemes based on dimensional power counting. 
Nevertheless, the essential difference is that non pertur- 
bative dimensional transmutation, i.e. the generation of 
dimension-full parameters not encoded in the potential, 
occurs due to the singular and attractive nature of long 
distance interactions already at the lowest order approx- 
imation consisting of OPE. This non-perturbative renor- 
malizability is the essential feature that makes this prob- 
lem particularly tough and so distinct from the previ- 
ous experience of perturbative renormalization on Effec- 
tive Field Theories or finite cut-off representations of the 
problem. 

The present work not only shows that the theoreti- 
cal requirement of renormalizability can be implemented 
as a matter of principle and as a practical way of con- 
trolling short distance ambiguities in the predictions of 
Chiral Perturbation Theory for the study of NN scatter- 
ing, but also that interesting physical and phenomenolog- 
ical insights are gathered from such an investigation. We 
have shown under what conditions such a program can 
successfully be carried out as a possible alternative and 
model independent way of describing the data by using 
very indirect, but essential, information on the implica- 
tions of chiral symmetry for the NN problem below the 
pion production threshold. 
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APPENDIX A: POTENTIALS 

For completeness we list here the potentials found in 
Ref. 3, and used in this paper. In coordinate space the 
general form of the potential is written as 

Vnn = Vcir) + fi ■ f2Wc{r) 

+ \ys{r) + Ti ■ f2Ws{r)\ a I ■ gi 

+ \yTir) + Tl • T2Wt('')] (3(7i • flT2 • f - (Tl ■ 172) 

+ [VLs(r)+n -faM^LsW] (Al) 

For states with good total angular momentum one ob- 
tains 

U]](r) = M {{Vc - iVs) + r{Wc - SW's)] , (A2) 
U]](t) = M{{Vc + Vs-yLs) 

+ t{Wc + Ws- Wls) + 2(^7 + rWT)\ , 

(A3) 

+ (j-1)(Vls + tW^ls) 
2j + 1 



(A4) 



6v^(jTT) 

2j + l 



M (Ft tWt) , (A5) 



2{j + 2 



2j 



^(Vt + tWt)] , 



(A6) 



with r = 2T(r + 1) — 3. Remember that Fermi-Dirac 
statistics requires (^—I'^^+s+T _ _^ 

The LO (OPE) potentials read {x — rriT^r ) 



y^OPE ^ 9^rn^ e " 
^ 487r/2 X 



(A7) 



all others being zero. 

The non-vanishing NNLO (TPE) potentials are given 

by 
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^ 3g'"^/ f 2ci + x2(l + + ^ + L + ^1 (6 + I2x + 10:r2 + Ax^ + x*) | 

^ ' 327r2/^ a;6 I V 16M ) ^ ^ ^ ^ 32M ^ \^ ^ -r ^^^^ y ^ ^ ^ ^ ^ M 

T^r-W = ^^{-(c, + ^)(l + .)(3 + 3. + .^) + j!^(36 + 72. + 52.^ + 17.3 + 2.^)}, 
^T^^^W = i28!r/4^4 { - 12i^o(2x) - (15 + 4x')K,{2x) + ^^^g^ (12^-^ + 24 + 20x + 9.^ + 2.^) } , 

4 5 

W^c''''(0 = [1 + 2.9'(5 + 2a;') - /(23 + Ux^)] Ki{2x) + x [l + lOg^ - /(23 + 4x2)] ^^(32;) , 

2 —2a: 

+ ^ ^'^^ — [2(3g2 - 2) (60;"^ + 12 + lOx + Ax^ + x^)] + g^x (2 + 4a; + 2x^ + ix^) | , 
Vr^'ir) = |!^{3xifo(2x) + (3 + 2x2)i^i(2x)-^^^^(6x-i + 12 + llx + 6x2 + 2x3)}, 



327r3/4 



p^TPiJ. ^ _ ff^Cg^ - l)m6 e-2^ , 



I + (1 + x)(3 + 3x + 2x^) - -|— (18 + 36x + 31^^ + Ux^ + 2x^) ] , 

L \ 4M / loAi J 



where iiTo and -f^i are modified Bessel functions. The 
NLO terms are obtained by dropping aU terms in 1/M 
and ci, C3 and C4. 



APPENDIX B: THE DIVERGENCE OF 
PERTURBATION THEORY FOR PERIPHERAL 
WAVES 

In this appendix we show that for a singular, attrac- 
tive or repulsive, potential at the origin which diverges 
like 1/r", there is always a finite order in perturbation 
theory where the phase shift diverges, regardless on the 
particular value of the angular momentum. Let us con- 
sider for simplicity the single channel case. The radial 
equation can be transformed into the integral equation 

poo 

ui{r) = jiikr) + GkAr,r')U{r')ui{r')dr' , (Bl) 



where Gk,i is the Green function given by 

kGk,iir,r') = Ukr)vi{kT')e{r' - r) 

+ Ukr')vi{kr)e{r -r'), (B2) 

where 9{x) is the Heavyside step function, 9{x) = 1 
for X > and 6{x) = for x < and j;(x) = xji{x) 
and yi{x) = xyi(x) are the regular and singular reduced 
spherical Bessel functions respectively. To regularize the 
lower limit of integration in Eq. (|B1|I one may assume 
a short distance regulator which will eventually be re- 
moved. The phase shift is given by 



iandi — — — 
k 



ji{kr)U {r)ui{r) 



(B3) 



In perturbation theory by successive iteration of Eq. IjBljl 
the Born series 



tan (5; = — - 



1 



dr 



l(kr) U{r) 



drdr']i{kr)U{r)U{r')GkAr,r')]i{kr') + 



(B 



is obtained. For our purposes of proving the divergence 
of perturbation theory it is sufficient to analyze the low 
energy limit. Using 5i — s- —aik'^''^^ and using known 
properties of the Bessel functions 



mix) 



(2^-1)!! 



(2;-f 1)!! 
The Green's function becomes 

- [21 + l)Go,i(r, r') = ^0(r' - r) + 



we get 



(B5) 



,1+1 



(B6) 



{2l + iyrai = / drr''+^U{r) 
Jo 

+ - -/ drr \ dr'{r'f^+^U{r)U{r') + 

21 -\- 1 In .In 



Since we only want to analyze the short distance be- 
haviour we can estimate the convergence of integrals 
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by using the finite range and singular potential U (r) = 
{R/r)^/R'^6{a — r). Thus, we see that in the first Born 
approximation the integral converges for 2/ + 1 > n — 2, 
whereas the second Born approximation requires 2Z + 1 > 
2{n — 2). This is obviously a more stringent condition. 
In general, at fc— th order convergence at the origin is 
determined by the integral 



finite order, even if the Born approximation was finite due 
to a high value of the angular momentum, I. 



APPENDIX C: LEADING SINGULARITIES IN 
THE SHORT DISTANCE EXPANSION 



dririU{ri) / dr2r2U{r2) . ■ 



The determination of the short distance behaviour 
(B8) from the full potentials is straightforward, but it is neces- 
sary to determine the number of independent parameters 
which is finite only for 2/ + 1 > k{n — 2), a condition in every channel and at any level of approximation. For a 
violated for sufficiently high k when n > 2. So, for n > 2 quick reference we list the leading singularity behaviour 
there will always occur a divergent contribution at a given in Table Hvl 
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TABLE IV: The leading short distance singularity of the NN reduced potentials, U = 2^iV , Eq. Q to LO, NLO and NNLO for 
all channels considered in this work. The sign of the coefficients for the one channel case (singlet and triplet) or the eigenvalues 
for the triplet coupled channel case determines the number of independent parameters. C3 = C3M and £4 = Mc4 are the 
dimensionlcss chiral constants. 
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